Effective Hamiltonians for periodically driven systems 
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The dynamics of classical and quantum systems which are driven by a high frequency (a>) field 
is investigated. For classical systems the motion is separated into a slow part and a fast part. 
The motion for the slow part is computed perturbatively in powers of to' 1 to order lu~ 4 and the 
corresponding time independent Hamiltonian is calculated. Such an effective Hamiltonian for the 
corresponding quantum problem is computed to order lu~ 4 in a high frequency expansion. Its 
spectrum is the quasienergy spectrum of the time dependent quantum system. The classical limit of 
this effective Hamiltonian is the classical effective time independent Hamiltonian. It is demonstrated 
that this effective Hamiltonian gives the exact quasienergies and quasienergy states of some simple 
examples as well as the lowest resonance of a non trivial model for an atom trap. The theory that is 
developed in the paper is useful for the analysis of atomic motion in atom traps of various shapes. 
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I. INTRODUCTION 

The interaction of cold atoms with strong electromag- 
netic fields results in many novel, interesting experimen- 
tal observations P, 0, ■ The relevant systems are char- 
acterized by an extremely high degree of control that en- 
ables one to explore various problems of general physical 
interest. The response to a rapid oscillating force is such 
issue and will be the subject of the present paper. 

Recently, in a series of experiments, atomic billiards 
were realized 0, • In these billiards atoms were con- 
fined by a standing wave of light to move in planes. The 
boundary of the billiards was generated by a laser beam, 
perpendicular to the plane of motion. This beam rapidly 
traverses a closed curve, which acts as the boundary of 
the billiard. The boundary of the billiard is assumed 
to be approximated by the time average of this beam, 
and the force applied by the boundary on the particles 
is approximately the mean force applied by the beam. 
One expects that this approximation is valid when the 
motion of the beam is fast relative to the typical veloci- 
ties of the atoms in the billiard. The billiards generated 
by the rapidly moving light beam motivated the present 
work. The more general physical problem, which is ex- 
plored here, is the description of classical and quantum 
dynamics in presence of fields that oscillate with high 
frequency 

In traditional atomic physics one typically assumes 
that the fields which affect the atoms have an amplitude 
which is constant in space and is time independent. This 
assumption is justified since the wavelength of the light 
field is much larger than the size of the atom and the 
electronic (internal) degrees of freedom react to the pe- 
riodic change in the field much faster then the external 
ones (center of mass coordinate and momentum). The 
main subject of traditional atomic physics is the response 
of the internal degrees of freedom to this field. Atomic 
spectroscopy is the most spectacular result of this line of 
research. The center of mass motion of the atom can be 



ignored in most laboratory experiments that explore the 
dynamics of the internal degrees of freedom. 

For the field of atom optics the effect of the internal 
degrees of freedom on the center of mass motion is im- 
portant, in particular near resonance of the external field 
with the internal motion (level spacing). The force on 
the center of mass due to the internal degrees of freedom 
is given approximately by a dipole force |2(. The sign of 
this force depends on the sign of the detuning of the light 
frequency from resonance (of the electronic levels). The 
motion of the atoms is manipulated by fields with am- 
plitudes which vary spatially resulting in a force on the 
center of mass of the atoms. In many cases the amplitude 
of the field can be assumed static. The atomic billiards 
described earlier consist of a time dependent field which 
results of the moving laser beam. Even at high frequen- 
cies of the motion of this beam one might expect that this 
time dependence will have some dynamical consequences. 
The question is most interesting when the wavelength of 
atoms is of the order of the size of the billiards. In this 
work the effect of a laser on the center of mass motion of 
the atoms will be modeled by a time dependent poten- 
tial. For some situations of physical interest this simpler 
model should still describe the dynamics in a high fre- 
quency field without the need to specify the dynamics of 
the internal degrees of freedom or the quantum aspects of 
the light field. Therefore, in the present work the atoms 
are modeled by point particles moving in a rapidly oscil- 
lating potential that varies in space. This description is 
relevant for a wide class of light-atom interactions and is 
not confined to models of billiards, which motivated the 
present work. 

The classical dynamics of particles influenced by a high 
frequency field was studied in several contexts. Kapitza 
investigated a classical pendulum with a periodically 
moving point of suspension [||. In this "Kapitza's pen- 
dulum" the motion can be separated into a slow part and 
a fast part which consists of a rapid motion around the 
slow part. The fast motion results in an effective poten- 
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tial for the slow motion. In some range of parameters this 
pendulum performs harmonic (slow) oscillations around 
the point where it points upwards. This point is unsta- 
ble in absence of the time dependent perturbation. This 
phenomenon is called "dynamical stabilization". Later, 
Landau and Lifshitz generalized this result for motion 
in presence of a rapid periodic force with a spatially de- 
pendent amplitude (see also and calculated the 
leading term in an expansion in powers of the inverse 
frequency. 

Dynamical stabilization is used to trap atoms in elec- 
tromagnetic fields. The most notable example is the Paul 
trap (3 • In this trap time dependent electric fields are 
used to localize ions in the region where the field am- 
plitude is minimal. The fields are well approximated by 
restoring forces which are linear in the distance from the 
equilibrium point. The resulting Hamiltonian is that of 
a time dependent oscillator. It is possible to find exact 
quantum mechanical solutions for this problem that are 
based on the corresponding classical system. That is, the 
states are simply related to the ones of the Harmonic os- 
cillator. Therefore the states of the motion in the Paul 
trap are known 0, . It is of interest to find some 
of the states of problems of a more general nature, even 
if only approximately. 

The work of Kapitza was first extended to quantum 
mechanical s yste ms in a pioneering paper by Grozdanov 
and Rakovic [Lj • They introduced a unitary gauge trans- 
formation resulting in an effective Hamiltonian that de- 
scribes the slow motion and demonstrated that its eigen- 
values are the quasienergies of the time dependent prob- 
lem. The effective Hamiltonian was calculated as an ex- 
pansion in powers of the inverse frequency. In that paper 
the analysis is restricted to a driving potential that has 
a particularly simple time dependence. Moreover, the 
final results are restricted to forces that are uniform in 
space, a situation natural in standard spectroscopy, but 
too restrictive for the interesting problems in atom op- 
tics. These restrictions are avoided in the present work. 

Other studies of quantum systems with periodic time 
dependent fields were also published. Gavrila 0, 
developed a perturbation theory for the Floquet states 
and the quasienergies in terms of the states of the time 
averaged problem. The scattering from a periodically 
driven barrier was studied by Vorobeichik et. al. |16| . 
by Bagwell and Lake 01 an d by Wagner [l^ , while the 
quantum and classical dynamics of some one dimensional 
systems were investigated by Henseler et. al. [lSj. In the 
limit of high frequencies the systems behave as if the 
particles were subject to an effective potential which is 
the time average of the time dependent one. Fredholm 
theory was used by Georgeot and Prange [2(J to study 
quasiclassical scattering from various systems, including 
a one dimensional periodically kicked potential. 

Another approach to time dependent systems (not nec- 
essarily periodic) is to use the Magnus expansion |21| in 
order to compute the propagator. Time periodic systems 
were used as examples in order to check the convergence 



of this expansion [22T I23I |2~I| . For these time periodic 
systems the Magnus expansion is of similar nature as the 
method presented here and the differences are discussed 
in Sec. UTTI 

There are numerous other works regarding periodically 
driven systems. Here we mention few of them. Of spe- 
cial physical interest is the ionization of atoms by light, 
(see |23 and references therein). Some toy models for 
ionization that consist of one dimensional time depen- 
dent S functions were treated rigorously [2(|. In par- 
ticular, it was shown that typically there is full ioniza- 
tion, namely, that at long times the probability to be 
at the bound state (of the time averaged problem) ap- 
proaches zero. For some arrangements of the S functions 
stable bound Floquet states exist. The transport through 
driven mesoscopic devices |27) and in the presence of os- 
cillating fields |23 attracted some interest. 

In the present work we study the dynamics of classical 
and quantum high-frequency driven systems. The clas- 
sical problem is discussed in Sec. [H] where the motion 
is separated into a "slow" and a "fast" part. A system- 
atic perturbation theory is developed for the motion of 
the "slow" part. The equation of motion of the "slow" 
dynamics is then computed to order w -4 , which is an ex- 
tension of the order uj~ 2 (resented in Ref. 0). This slow 
motion is shown to result from an effective Hamiltonian. 
In Sec. lIIII an adaption of the Floquet theory to the prob- 
lem is reviewed. An effective (time independent) Hamil- 
tonian operator is defined following and generalizing |13| . 
The eigenvalues of this operator are the quasienergies of 
the system. This effective Hamiltonian is then computed 
perturbatively (to order lu~ 4 ) in Sec. II VI The restrictions 
introduced in |l3j are avoided and consequently detailed 
expressions for the various terms of the effective Hamil- 
tonian are calculated explicitly. The classical effective 
Hamiltonian of Sec.[H]is found to be the classical limit of 
this quantum effective Hamiltonian. Two known exactly 
solvable simple examples of the method are presented in 
Sees. Ivl and I VII In Sec. IVIll the scattering from a time 
dependent potential is discussed. In particular, the reso- 
nances of the time dependent problem are found to agree 
with those of the effective (time independent) Hamilto- 
nian of Sec. II VI Finally, the results, the implications and 
some related open problems are discussed in Sec. IV111I 



II. CLASSICAL MOTION IN A HIGH 
FREQUENCY POTENTIAL 

In this section the dynamics of a classical particle mov- 
ing in one dimension under the influence of a force which 
is periodic in time is studied. Typically, solutions for 
time dependent problems can only be attained numer- 
ically. However, when the period of the force is small 
compared with other time scales of the problem it is pos- 
sible to separate the motion of the particle into "slow" 
and "fast" parts. This simplification is due to the fact 
that the particle does not have the time to react to the 
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periodic force before this force changes its sign, namely, 
the contribution of the periodic force to the acceleration 
in one period is negligible (compared to the contribution 
of the effective force in a sense that will be specified in 
what follows). Thus we will consider the limit of small 
periods (or large frequencies) of the driving field. 

The leading order (with respect to l/u>) of the dy- 
namics was computed by Kapitza 6] for the "Kapitza's 
pendulum" , namely a pendulum where the point of sus- 
pension is moved periodically. It turns out to be very 
general Q. Here the next order is computed and it is 
demonstrated that the equation of motion of the slow 
part of the dynamics can be derived from a time indepen- 
dent Hamiltonian. This Hamiltonian will be computed 
explicitly to order l/w 4 . Later, this Hamiltonian will be 
compared with an effective Hamiltonian which will be 
derived for the corresponding quantum problem. 

The existence of such a Hamiltonian might seem to 
contradict the fact that the time dependent dynamics do 
not possess a constant of motion. Moreover, the classical 
motion may be chaotic. The existence of this effective 
time independent Hamiltonian implies that a constant of 
motion exists for the slow dynamics (it is just the effec- 
tive Hamiltonian) , and for a one dimensional system the 
slow dynamics is integrable. To avoid confusion, it should 
be emphasized that the effective Hamiltonian depends on 
a coordinate which describes the "slow" part of the mo- 
tion. This coordinate is not the location of the particle 
(although they are almost identical at high frequencies). 
The actual motion consists of rapid motion in the proxim- 
ity of the trajectory of the slow dynamics. The relation 
between the slow coordinate and the coordinate of the 
particle is nonlinear and extremely complicated as will 
be demonstrated in what follows. We will demonstrate 
that an effective Hamiltonian for the "slow" motion may 
exist. 

Newton's equation for the motion in the periodic field 
is given by 



where V\ is a periodic function of tot of period 2-k and 
its average over a period vanishes. We denote derivatives 
with respect to coordinates by primes and with respect 
to time by dots. This separation of the potential to an 
average part, Vq(x), and a periodic part with vanishing 
average, Vi(x, tot), is natural and will simplify the follow- 
ing calculations. We look for a solution of the form 



x(t) =X(t)+£(X,X,ut) 



(2) 



where X — ^Sr and 

dt 



2tt 



£=— / dr^X,X,r)=0. (3) 

Z7T 



The bar denotes in this paper the time average over one 
period. The fast part of the motion, that is nearly peri- 
odic in time, is denoted by £. It will be shown later that 
it can be chosen to depend only on X and X, but not on 
higher order time derivatives. Since X and X are slowly 
varying functions of the time t, £ is not periodic in the 
time t, in spite of @. The coordinate X describes the 
slow part of the motion and its equation of motion will 
be computed in the following. Our method of solution 
is to choose £ so that Q will lead to an equation for X 
which is time independent. An exact solution using Q is 
too complicated to obtain. However, at high frequencies, 
one can determine £ order by order in 1/uj. In order to 
separate terms in powers of the frequency it is convenient 
to introduce the new time variable r = ut. Using 



Tt- UJ d^ + dx x + ^x X 



(4) 



dt 2 



■Vi(x)-V{(x,ut), 



(1) Newton's equation is given by 



X 



dr 2 



2lo 



d 2 z 

dXdr 



X 



d 2 (j 

dXdr 



X 



dX 



X 



dX 



X 



d 2 z 



X" 



d 2 i ■ •• 

2 ^11 



d 2 i 



x 2 



dX 2 ' ' dX8X" ~ ' dX 2 . 
= -V{(X + Z)-V{(X + £,t) 



(5) 



The variables r and t will be treated as independent vari- 
ables. This calculation is similar to the ones performed 
within the method of multiple time scales. Indeed, the 
result of the following calculation is equivalent to the 
one obtained using the method of multiple time scales, 
as demonstrated in App. [S] 

In the limit of high frequencies £ is going to be small 
(of order tu~ 2 ) and therefore it is convenient to expand 



Vo[X+£) and Vi(X+£, r) in powers of £ (we assume that 
Vq and V\ are smooth functions the coordinate). Then £ 
is expanded in powers of 1/uj 

oo 1 
i=l 

The £j are chosen so that the equation for X that results 
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from (0 does not depend on r. 

Before obtaining the slow equation of motion from J3J , 
order by order, there are two points regarding our method 
of solution which should be discussed. First, we note that 
the fast part £ is expanded in powers of 1 /to while X is not 
expanded, which seems to be inconsistent. One may also 
expand X in powers of 1/ui as X — J2°l -\Xi. When 
one does so the equation of motion for X is replaced by 
a series of equations for Xi. In this series of equations 
each Xi can be determined from the lower order terms 
Xj, where j < i. This is the standard method of sep- 
aration of time scales and its application to the present 
problem is demonstrated in App.[S] These equations are 
equivalent, in any order, to the equation of motion of the 
(unexpanded) X which will be obtained in what follows. 
At a given order co~ n of the present calculation all con- 
tributions that are found by the method of separation 
of time scales are included, but some of the higher order 
terms are included as well. Second, we note that while we 
assumed that £ depends only on X and X, higher order 
derivatives of X with respect to time appear in equation 
(jSJ). In the leading order in 1/w, as will be demonstrated, 
one can replace X by — —Vq(X). The error is of higher 
order in l/w, leading to the correct contribution to 
at the order where X appeared. Corrections of higher 
orders of 1 juj to £i result from the corrections of higher 
orders to X. These corrections will affect the £, with 
j > i, since these are chosen to cancel the r dependence 
at any given order. Higher order derivatives of X can be 
found by repeated differentiation of X. This enables us 
to obtain an expression for £ that depends on X and X 
but not on higher order derivatives of X. 

To proceed we gather all the terms in equation 
using (JHJ, that are of the same order, say ui~ n , and choose 
Cn+2 (which is still undetermined), so that the explicit r 
dependence cancels. In the leading order, (w), the only 
contribution is 



9 2 6 



dr 2 



0. 



(7) 



Therefore we can choose 



a = o. (8) 

In the next order (w°), we find the contributions 



X 



d^2 

dr 2 



X(X)-V{(X,t). (9) 



Our goal is to balance the r dependence. To do this we 
have to solve 



8t 2 m 



(10) 



moreover, we also require that £2 is periodic in r. The in- 
tegral over the right hand side (RHS) of equation ljTU|) can 



have terms which are time independent and thus £2 can 
grow linearly in r. To ensure that £2 is small even at long 
times such secular terms must be avoided. This can be 
done by requiring that the time integral has a vanishing 
average over a period. Let f(x,r) be any periodic func- 
tion of r with a vanishing average, / = 0. Assume that 
the Fourier expansion of / is given by / = fn emT > 

then we define the following integral, 



f [/] = ]T -f n e 



(11) 



n#0 

and its repeated application will be denoted by 

c(2)t 



(12) 



and j applications by 



[/] = / 



[/]• 



(13) 



j times 



This definition, which is actually a specific choice of the 
integration constant, is natural since it ensures that the 
result is periodic even after repeated integrations. It also 
helps to separate periodic terms (with vanishing average) 
and secular terms (which will be time independent in the 
current calculation). Integration of l(TT)|) implies 



1 r i2)r 

6 = --/ [V((X,r)] 



(14) 



Note that we did not really find the general solution 
of ©, but rather chose £ so that it is satisfied. Substitut- 
ing £2 in equation gives the leading order equation for 
the slow coordinate X. The terms in this equation are 
just the time independent terms which were not canceled 
by 6, 



mX = -Vq(X). 



(15) 



The contributions from JSJ at the next order, oj , are 



n 3 . ~v, a 2 

<9r 2 



2X 



dXdr 



(16) 



which, with the help of 114|l . is satisfied by 



o /-(3)t 

& = / [V{'{X,t)\. 
m J 



(17) 
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The terms of order w 2 in equation (J3J) are given by 



d 2 U 
dr 2 



2 X 



nXih ! " i)Xdr) 



a 2 6 



+ x 



X 



56 
dX 



X 



dX 2 



X'(x)^-v;'(x,r^ 2 . 



(18) 



Substituting and JJZJ leads to 



^4 
9r 2 



2 



(2)T 



T/" 



m- 



(2) 



[y/] - -x' 

171 



(2) 



Fi (3) (*,r) 



3X 



m 



(2) 



(19) 



Equation (|19|l cannot be solved, if £4 is required to be 
periodic in r, since the RHS has a non vanishing average 
which will lead to solutions that grow like r 2 (these are 
the secular solutions that one wishes to avoid when using 
multiple time scales analysis). We will choose £4 so that 
it will balance the r dependent part of l|19|) and will be 
periodic in r. The remaining r independent terms in 
(|19f) will be included in the equation of motion of the 
slow coordinate X. Defining 



fi(X,r) 



1 r(2)r -I A2)t 



(20) 



order w~ 5 . To obtain the next correction to equation H22|) 
one needs only the average over r of the terms of order 
lu~ 3 . The reason is that £5 will be chosen in such a way 
that it will cancel all the periodic terms with vanishing 
average. This further simplifies the calculation since all 
the terms (except ml) on the LHS of JSJ) have vanishing 
average (over r) , thus only terms from the RHS can 
contribute to the equation of the slow coordinate. In 
this order the contributions to the equation of the slow 
coordinate X can result only from 



TO ~ TO 



(23) 



and choosing 

£4 = 4 / [v{] 



m 

nv /-(4)r 



(2)r 



LA 



3X 2 



ir2 r (4)T 



(3) 



(21) 



balances all the r dependent terms on the RHS of i|19|) 

but leaves an extra term 



1 fW T 

-v" 

m J 



[V{] 



which is not balanced. This is actually a term of order 
that is left on the RHS of (0 when we substitute 



£ in ©. 

then 



The resulting equation for the slow motion is 



mX = -Vq{X) 



moj 



-VI 1 
2 v i 



(2)r 



[Vn + Oiuj- 3 ). (22) 



This is the leading order correction due to the periodic 
potential V\. It was calculated before H0. With the 
help of (|15(l or of the leading term in (|22|l X can be 
eliminated from the expression d2U for £4. This method 
allows us to compute corrections order by order. We will 
continue the calculation up to order cj~ 4 . 

The next order is ui~ 3 . We do not need to compute £5 
explicitly since it can only change the slow equation in 



The first term will vanish since Vq is r independent and 
£3 has a vanishing average. The second term can be 
computed using (fTTIl 



2X 



V>> 



(3)r 



[V{'] 



[V'»]=0. (24) 



In the last calculation we have used integration by parts 
and then the fact that the average of a derivative of a 
periodic function over a period must vanish. This leads 
to the conclusion that one can choose a periodic £5 in 
such a way that all r dependent terms of order o>~ 3 in 
© are canceled. 

We turn to the order uj~ 4 that is the last order that will 
be considered here. Again one can get the contributions 
to the equation of X by averaging terms of this order 
in JSJ). The average over r of the LHS vanishes and the 
contribution of the terms on the RHS is 



TO 



1 



n v S2 



TO 



-V (S) £ 2 
2 1 42 ' 



(25) 



The first term will vanish but the other terms have a non 
vanishing average. Using (|14|) . (|21|l and integration by 
parts (in the averages) yields 
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In the last term the X in £4 was replaced by — Vg/m resulting in errors that are of order cu~ 6 in the final result. Eq. 
(|26[1 gives the o> -4 contribution to the equation for the slow coordinate X. 

The equation for X to order w -4 is obtained when £ is substituted into (JSJ and the remaining terms are averaged 
over r resulting in 



1 ~r r 1 n\( r (2)T V 1 r {2]T r {2)T 

mx = - — 2 J in] J [vn - { J raj - wn J [vn 

1 mf f (2)T Y 1 7^ 7^ T iX 2 !^ f (2)T r ^1 

-2^ () (/ K'i)-^ y '7 K1 / [V?, W ray K1 

- S(/ W >,<<])W»). 

It is instructive to introduce the effective potential 



V eff (X) = V + 



1 



2muj 2 



[vn 



r(2)T 



VT 

2mV 1 



[^] 



»(2)r 



2m 2 i 



— V" 

,4 ^0 



(28) 



Substituting l|28|) in (|27|l results in the equation of the 
slow motion 



mX = — 




Given a solution of this X(t), the solution for the original 
problem can be easily obtained (to the appropriate order 
of 1/ui) since £ is known in terms of X (see lfH)l. l(T7j) 
and H21|l). From these equations one sees that in the case 
where the oscillating force V{ is independent of position 
X the fast coordinate £ is independent of X and X to 
the order oj~ 4 (note that in H21(l only the order ui° of 
X is required see also JA161 ). The final result of is 
confined to the case where V{ is independent of X. The 
equation l|29|) can be derived from the Hamiltonian 



P 2 3 / r(a)T y 

ft W - si+^w+^y raj 

+ 0(^- 5 ) (30) 
where P is the momentum conjugate to X . 



We have shown that using the natural separation of 
time scales it is possible to separate the motion of a par- 
ticle in a high frequency periodic field into "slow" and 
"fast" parts. The slow dynamics can be derived from 
an effective Hamiltonian which is time independent. We 
turn to discuss the corresponding quantum problem. 

III. FLOQUET THEORY AND THE EFFECTIVE 
HAMILTONIAN 

Consider a quantum system with a Hamiltonian that 
is periodic in time, H(t + T) — H(t). Such systems 
can be treated using Floquet theory jH M, HI M, HI 
The symmetry with respect to discrete time translations 
implies that the solutions of the Schrodinger equation 

ih-^iP = H^, (31) 

are linear combinations of functions of the form 

ipx = e~ l ^u x (x,Lut) (32) 

where the u\ are periodic with respect to ut with period 
27T, that is u\(x, u>{t + T)) = u\{x,ojt) with to = 2ir/T. 
The states u\ are called the quasienergy or the Floquet 
states and A is referred to as the quasienergy (we will also 
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call the states ip\ quasienergy states). This is the content 
of the Bloch-Floquet theorem in time. The states u\ are 
the eigenstates of the Floquet Hamiltonian 



this equation is 



d 

-ih- + H. 
at 



(33) 



The quasienergy (or Floquet) states have a natural sep- 
aration into a "slow" part e~ l ^~ (with the natural choice 
< X/h < lu), which includes the information about 
the quasienergies, and to a fast part u\(x,uit) that de- 
pends only on the "fast" time u>t. It is expected that one 
will be able to find an equation of motion for the slow 
part of the dynamics as was done for classical systems in 
Sec.[n] Such an equation will include information regard- 
ing the quasienergies of the quantum system, and will be 
developed in what follows. It establishes a natural link 
between the separation of the fast and the slow motion 
in classical mechanics, which can be formalized by the 
theory of separation of time scales, and Bloch-Floquet 
theory in quantum mechanics. 

It is known that one may write the propagator in the 
form 33] 



U(t) = P(t)e 



-itg/n 



(34) 



where Q is self-adjoint and P is unitary and periodic with 
the period of the Hamiltonian. The eigenvalues of Q are 
the quasienergies of the system provided that the eigen- 
states of Q are in the domain of H(0). Sometimes Q is 
called the quasienergy or Floquet operator. The actual 
calculation of Q might be complicated. Such an operator 
was calculated in [22j and in [13| by introducing expan- 
sions for P and Q. The result turns out to depend on the 
phase of the periodic part of the Hamiltonian or on the 
initial time. (See for example [22j, equations (25) and 
(26) and [3, equation (16)). Inspired by an 
approach of a somewhat similar spirit is used. 

The goal is to find a unitary gauge transformation 
e* F M, where F{t) is a hermitian operator (function of 
x and p) defined at a certain time t, which is a periodic 
function of time with the same period as H, such that in 
the new gauge the Hamiltonian in the Schrodinger equa- 
tion is time independent. Such a Hamiltonian was found 
by Grozdanov and Rakovic |13| if the time dependent 
part is of the restricted form Vqr = V{x) sm{uit + 9). 
It was analyzed with the further strong restriction, that 
for one dimension takes the form W- — const (uniform 
force). In what follows a general analysis that is free 
of these restrictions is presented. Applying e lF to both 
sides of equation (|31|l and adding ifi (j^e 1 ^ i> to both 
sides leads to 

ih^ (e^v) = e iP Hff> + ih (§-/ P ) i>- (35) 



ih di ^ = G ^ 



where the Hamiltonian is 



G = e lF He- tF + ih 



de l 



Of 



In the classical limit it reduces to 

OF 



G = H — h 



dt 



(36) 



(37) 



(38) 



Therefore in the classical limit — TiF is the generating 
function of the canonical transformation corresponding 
to the unitary transformation e~ lF [34|. 

Let us assume that such an operator F exists so that 
G is time independent. Then the cigenfunctions of G are 
v\(x), their evolution takes the form 

4> x {t,x) = e- i %v x {x). (39) 

These states, in the original gauge, correspond to 

tl>x{t,x) = e~ iF 4> x = e-^e- iF v x {x), (40) 

since F does not include any time derivative. The func- 
tion e~ lF v\ is periodic in time with the period of H and 
therefore ip\ is a Floquet state with quasienergy A (mod 
flu). It should be compared with (|32l) with the identifi- 



cation U\ 



-iF 



v\. It is assumed that e (and e ) 



are such that they map the domain of H(t) into that of G 
and vice versa. This may not be true in general, and one 
cannot exclude the possibility that examples, where only 
some of the quasienergies can be found using this method, 
exist. For example, problems of this nature may occur if 
for a function ip in the Hilbert space of H , the function 
e lFr ip is not in this space. The limitations on the validity 
of the method should be subject to further mathematical 
studies. 

To emphasize the difference between the effective 
Hamiltonian G and Q given by equation (|34l) let us write 
the propagator in terms of F and G. To propagate any 
state in time using G it has to be transformed to the 
time independent gauge, then propagated and then trans- 
formed back. This results in the propagator 



U(t) 



= p -iF(t) p -i^ p iF(0) 



(41) 



In terms of the functions in the new gauge, 



Since G and F do not commute, G generally differs from 

a of (J3H). 

We note that an approximate solution of the time de- 
pendent problem in terms of an expansion of F and G 
has some superior properties compared to the more cus- 
tomary expansion of V and Q of (|34() . For instance, if F 
is hermitian at any order than e %F is manifestly unitary 
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while some care is needed to obtain unitary approxima- 
tions for V . In addition, G does not depend on the phase 
of the time dependent field while Q does depend on this 
phase (see [11122^. Therefore in the present work a de- 
scription in terms of G and F is used rather than one in 
terms of V and Q. 

In the following section the derivation of G and F will 
be presented explicitly as an expansion in powers of 1/ui. 
It will be shown that at high frequencies F can be cho- 
sen to be small, of the order of 1/uj. In this limit one can 
easily calculate matrix elements of an observable O be- 
tween quasienergy (Floquet) states using the eigenvalues 
and eigenstatcs of the effective Hamiltonian G 



i , 



2 W 



A 2 ) +«(0A 

F, 6 



F,0 



F 



(42) 



The result is an effective expansion in powers of 1/lu. It 
was obtained with the help of l|B3(l presented in App. [5] 
Since observables have a meaningful classical limit h — > 
their expectation should reduce to the expansion in pow- 
ers of 1/lu for the corresponding classical quantity as 
calculated in Sec. [H] and App. EJ The expansion of G 
presented in the next section can be considered an ex- 
tension of the multiple scales analysis to quantum me- 
chanics. The effective Hamiltonian that will be obtained 
will be compared with the classical Hamiltonian for the 
slow motion that was computed in Sec. ^ 



IV. THE EFFECTIVE HAMILTONIAN OF 
QUANTUM SYSTEMS WITH A HIGH 
FREQUENCY POTENTIAL 



In Sec. IIIII we demonstrated that the quasienergies and 
Floquet states of a quantum system can be determined if 
one can find a gauge transformation so that the Hamil- 
tonian is time independent. The transformation and the 
resulting effective Hamiltonian are obtained here. Typ- 
ically F and G cannot be computed exactly. For high 
frequencies one can determine F and G order by order in 
1/lu. In the following we present a derivation of F and 
G accurate to order 1/aA 

We consider the Hamiltonian (that is more general 
than the one studied in ^|) 



H 



P 

2m 



+ V Q (x) + V 1 {x,Lut). 



(43) 



This is the quantum system which corresponds to the 
classical system that was discussed in Sec. [H] It should 
be noted that the method which is described in the 
present section also applies to Hamiltonians that differ 
from l|43|) , for example in presence of magnetic fields and 
for spins (see Sec. I VI \ . We choose to examine the Hamil- 
tonian (|43[1 since it is of interest to compare the resulting 



effective Hamiltonian with its classical counterpart 13011 . 
As mentioned in Sec. IIIII we are looking for a unitary 
transformation e lF so that the resulting Hamiltonian l|37|l 
is time independent. It is convenient to define r = tut, 
since the Hamiltonian depends on time only through r. 
Using this definition (|37|l is given by 



G 



He 



-iF 



■ iTiuj 



de 



iF 



8t 



(44) 



At high frequencies F is assumed to be small, of the or- 
der of 1/lu. An assumption that will be explicitly satisfied 
by the following calculation. This enables us to expand 
G and F in powers of 1/lu and to choose F so that G 
is time independent in any given order. The expansions 
are given by 



g = y —d n 



n=0 



and 



oo 



(45) 



(46) 



The periodicity F(t + 2ir) — F(t) is assumed. The 
calculation is performed by computing G; in terms of 
Fx, ■ ■ ■ , Fi + i and then choosing Fi + i so that Gi is time 
independent. The terms in (|44|l are calculated with the 
help of the operator expansions (presented in App. !b|) . 



e iF He -iF = H + i 



I 

3! 



F,H 
F 



1 

2! 



F. 



F, H 



follows directly from i|B3|) and 



de 



iF 



-iF 



Or 



OF 



F, H 



1 

2! 



(47) 



F, 



OF 



i 

3! 





' df 




F, 







(48) 



is a result of (|B5p and (|B6ll . 

In the leading order, O(oj°), Gq is given by 



V 2 - - dF\ 

Go = j- + V (x) + V^x, t) - h—±. 
2m or 



(49) 



The potentials Vq and V\ do not depend on p. To cancel 
any time dependence we choose 



Vx{x,r) 



(50) 



It is easily computed in the coordinate representation. 
Note that F\ is determined only up to a hermitian time 
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independent operator. It was assumed to vanish here. 
Substituting in gSJl leads to 



Go = f - + V (x). 
Im 



(51) 



This is the leading order of the effective Hamiltonian. 
The dynamics do not depend on the fast time depen- 
dent potential V\ as expected. Corrections due to V\ will 
appear at higher orders in 

At order 1/lu the effective Hamiltonian obtained from 

GU-EHJ is 



Gi 



Fi,H 



dF 2 
l 97 



ift 
~2 



Fx 



dFi 



(52) 



Note that Fi , given by (|50|) , depends only on the coordi- 
nate and therefore it commutes with its time derivative 
and also with Vq. If a periodic F 2 can be chosen so that 



m 

dr 



Fx,H 



Fx 



P 

2m 



then G± vanishes. Indeed, by choosing 



i f (2)T i 
2m / m 



(2)r 



we obtain 



d =0. 



(53) 



(54) 



(55) 



We have presented F 2 in the coordinate representation 
since the simple dependence of H on momentum makes 
it the most natural representation. We will use it also 
when calculating higher orders. 

At the next order, uj~ 2 , G 2 found from is 



G, = 



F 2 ,H 



Fl 



Fi,H 



dF 3 



itl 

~2 

h 
6 



- dF 2 



iTi 
~2 



F 2 , 



dF\ 
dr 



Fi 



Fi, 



dF\ 
dr 



(56) 



Substituting H = Go + fo^pjr and using (|53[1 to eliminate 



the commutation relation 



F U H 



results in 



G, = 



F 2 , Gc 



- h 



dF 3 
dr 



iTi 
~2 



F 2 



dFi 



(57) 



We can choose a periodic F 3 in order to balance the time 
dependence of G 2 . Note that G 2 has some time inde- 
pendent part that cannot be canceled by a periodic F 3 . 
Therefore we separate G 2 into a r independent part and 
a part that is periodic with vanishing average and choose 



F3 so that the latter vanishes (in 15711 ). For this purpose 
F3 must satisfy 







'fi 




fi 9Fi~ 




F 2 , Go 


< 






F2 ' ~d7 


) 



dF 3 _ i 
~ h 

'(58) 

where Go is given by l|51|) and an average over a period is 
denoted by bar. After some algebraic manipulations F 3 
is found to be 



F 3 = 

) - 



h 

rrv 

h 

Am 2 
1 

2mh 



(3)r Q2 n 

a_ n_ 

11 J ™2 



(3)r 



(3)r 



V, 



(A) 



dx 2 m 2 



(3) 



d_ 

O.r 



[^i] + / 3 (£,p) 



(59) 



where 



imh 



dF x 



F 2 



"a7 



,(2)t ,(2)t 



(60) 



The constant of the integration over r is the hermitian 
operator f 3 that depends only on x and p, and will be de- 
termined at the next order. Later we will use the freedom 
to choose f 3 to cause G 3 to have simple form. Using Ij59(l 
in H57[) will cancel the time dependent terms resulting in: 



G, = — 



F 2 , 



~d7 



1 , [ {2)T 



[V{] 



(61) 



where we have used integration by parts. 

The calculation of G 3 and G4 can be performed along 
similar lines. Since it is tedious it will be outlined in 
App. [C] It starts from (|Clfl and (|C9|1 that are obtained 
from l|44|) - (|48|l . Using the freedom in the choice of f 3 we 
choose it to satisfy (|C4|) . so that G 3 vanishes. Then F4 is 
found to satisfy (|C6(1 , that is integrated to take the form 
(|C7|) . The time independent part of F4 is denoted by 
Using the freedom of the choice of gauge, we choose f± 
so that in the classical limit the effective Hamiltonian G 
reduces to its classical counterpart (|30|) . This results in 
G4 of i|C22|) . In App. |0 we calculated also F4 that was 
required for the calculation of G4 but not F5 required for 
the calculation of G5. 

The freedom in the choice of gauge was used here and 
the time independent parts of the Fi were chosen in a 
specific way. Generally this choice is arbitrary. In the 
present work, a choice was made so that in the classi- 
cal limit the effective Hamiltonian reduces to the specific 
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classical counterpart 1|30|) , that resulted in a natural way 
within the derivation of Sec. |H] In 13], on the other 
hand, the choice of the time independent parts of the Fi 
is made so that the average of the fast variables over a 
period reduces to the slow variables within an order uj~ a 
calculation. It is found there, that with this choice, the 
requirement can be satisfied only if the oscillating force 
is independent of position (in one dimension). 

We have used perturbation theory (in l/u>) to ob- 



tain a periodic gauge transformation e lF and an effective 
Hamiltonian G so that the quasienergies are the eigenval- 
ues of G. Its eigenstates are related to the quasienergy 
states by Q4U[1. For a Hamiltonian of the form l|43|) this ef- 
fective Hamiltonian is given by equations l|44|l . (|45|) . I|51|l . 
iJSSJl, lEU), (|U5|l and l|C22jl . Collecting all contributions 
one finds, 



d = hi + fe// + 4^ (P 2 ^ + 2 P9{x)p + g(x)f) + ^V q + 0(u;- 5 ) 



where 



V Bff (x) = V Q (x) + 




2mu> 2 



1 / 

Vf / [V{] 



2m 2 w 4 



(62) 



(63) 



is the effective potential corresponding to l(2"%|) , 
3 



2m 3 




[V{> 



(64) 



is the coefficient of P 2 in H3()(l while 
1 



K2)r r 



v„ 



8m 3 



(•'!) 



(65) 



is a quantum correction to the classical Hamiltonian (its 
form obviously depends on the ordering of operators in 

153). 

The effective Hamiltonian is the main result of this 
section. The classical limit of l|62[) is the classical effective 
Hamiltonian (|30() . The freedom of gauge in the quantum 
problem was used, and fa and fa were chosen specifically 
to achieve this. We did not use the freedom of a canonical 
transformation in the classical calculation. The specific 
canonical transformation from the Hamiltonian (|43ll to 
the Hamiltonian (|30|l is generated by the classical limit 
of — fiF with the specific choice of the time independent 
parts that was made in the present work. 

The perturbation theory that was developed here en- 
ables one to calculate not only the quasienergies but also 
the corresponding Floquet states. If the eigenfunctions 
of G are known, then the quasienergy (or Floquet) states 
can be computed up to order u;^ 4 using equation 14011 
with 



i 



(2)7 



where F3 is given by (|59|l and i|C4(l while F4 is given by 
EJ , iHTQll and Oil . 



V. THE HARMONIC OSCILLATOR DRIVEN 
BY A PERIODIC EXTERNAL FORCE 

In this section a simple example that may help to clar- 
ify the meaning of the effective potential and the gauge 
transformations used in Sec. IIVI is discussed. It is the 
harmonic oscillator driven by a force that is periodic in 
time. It is defined by the Hamiltonian (|43|) with 



and 



V (x) 



Vi(x, 



1 



2 2 
muj n x 



£xcos(t). 



(67) 



where m is the mass of the particle, lo q is the classical 
frequency of the oscillator and £ is the amplitude of the 
driving force. We assume the non resonant case uj ^ 
loq. For this simple system one is able to compute the 
quasienergies and the Floquet states exactly. These were 
computed in several previous works 0, 0> 13 HE HE 

m. 

The high-frequency perturbation theory of Sec. IIVI can 
be used in order to calculate F and G. The calculation 
is straightforward and substitution of (|43|l with (|67ll in 
(E2J> with (ESJl-lESIl results in 
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and substitution in Il66[l yields 



£x 

F = — sm r 
nuj 



£% 



~~2 C0ST 7T 
r ox 



8mhuj 3 



sin(2r) 



ha; 3 
d 



SUIT 



—^-j cos T — 

mui^ ox 



+ O (uj- 5 ) (69) 



In this case both G and F have a simple form. The effec- 
tive Hamiltonian G reduces to a Hamiltonian of a time 
independent Harmonic oscillator. Equation I69|) implies 
that the gauge transformation is of the linear form (to 
order 1/w 4 ) 



F = A(t)x + B(t)+C(t)p 



(70) 



where A, B, C are real periodic functions of time. It turns 
out that l|70[) is exact with 



A(t) 
B(t) 



h(u> 2 - uJq) 



— j-= sin(2wi) 

£ 

C W = r i 2 2TCOs(wt). 



(71) 



For this purpose one calculates and the commutators 
required for the expansions (|47() and 148(1 . These are very 



simple since only 



F,H 



and 



F,H 



do not vanish. 



Requiring that the linear terms in x and p in the effective 
Hamiltonian G vanish determines A(t) and C(t) of (|71ll . 
and the requirement that it is time independent results in 
B(t) of <|71jl . Finally the Hamiltonian G takes the form 



G 



2m 



2„2 



£ 2 



4m(w 2 — u/q) 



(72) 



Expansion of (|72|l in powers of 1 juj results in l|68[) while 
the expansion of (|70|l with l|71|l results in (|69|) . The eigen- 
values of G are the quasienergies of l|43|> with l|67|) . Since 
G is a harmonic oscillator Hamiltonian they are given by 



(73) 



with non negative integer n. Note that these quasiener- 
gies are determined only up to integer multiples of the 
period of H, that is, they are given modulo Ticu. The 
Floquet states of this oscillator can also be computed. If 
v n (x) is the textbook eigenstate of the harmonic oscil- 
lator l|72(l . with the energy E n , than the corresponding 

Floquet state is ipnix^) = e~ l ~%~ e~ % ^ v n (x). 

The Floquet states (u\ of l)32|l ') are defined only up 
to powers of e muJt , that is, multiplying one of them by 
such a factor will also give a quasienergy state with the 
quasienergy shifted by nhu. It is possible to verify that 
these states are indeed eigenstates of the Floquet Hamil- 
tonian , and their eigenvalues are given by (|73|) . One 



may require that all the quasienergies are in the interval 
[0, hu>) and describe the dynamics by such states. If uj/ujq 
is irrational almost any value in this interval corresponds 
to a quasienergy. If uj/luq is rational there is only a finite 
number of (infinitely degenerate) quasienergies. 

This harmonic oscillator driven by a periodic force 
serves as a simple example for the method developed in 
Sees. HTll and llVl For this simple example one can com- 
pute the quasiener gies and Floquet states of the system 
exactly [ill 111 lillBl IM l37| . The perturbation theory 
described in Sec. |^ leads to an expansion in powers of 
1/uj. So far, we did not discuss the convergence of the 
series for the effective Hamiltonian G and for F. For the 
system discussed here both series converge when to > luq . 
The series can be resummed also for w < u>o and the re- 
sult is valid for all w ^ loq. In general the convergence 
properties of the series for F and G are not known. It is 
possible that the results of Sec. II VI are meaningful also 
for some cases where the series do not converge, if these 
can be resummed, like the ones of the present example. 



VI. DYNAMICS OF SPINS IN TIME 
DEPENDENT MAGNETIC FIELDS 

Another simple example that demonstrates the meth- 
ods presented in this work is a spin in a field, which 
is a combination of a static and a periodic time depen- 
dent magnetic field. This example demonstrates that also 
Hamiltonians that are not of the form p 2 /2m + V(x) can 
be treated in the way presented in Sees. Hill and IIVI The 
systems that are considered here consist of a spin in a 
constant magnetic field combined with a perpendicular 
periodic field with linear or with circular polarization. 
The Hamiltonian for the linearly polarized field is given 
by 



Hi = -ojqI z + lui cos(u;t)I x 
while for the circularly polarized field it is 



(74) 



-ujqIz + lui (^cos(ujt)i x + sin(cijt)ij / ^ . (75) 



For a spin in a circularl y p olarized field the problem was 
solved exactly by Rabi [33. This system also appears in 
textbooks as a paradigm of time dependent two level sys- 
tems j3^ . Our goal is to demonstrate that for this simple 
system, the quasienergies can be computed exactly using 
the method presented in Secs. lIIII and llVI First we derive 
some results that are valid for any Hamiltonian linear in 
the spin operators. 

These spin problems turn out to be simple since the 
spin operators have a closed algebra, 



1% 5 



ihl. 



ihL r 



ihL. 



(76) 

The effective Hamiltonian (|44|) is obtained with the help 
of the expansions in commutation relations (|47|l and 14811 . 
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For a Hamiltonian and F that are linear in the spin op- 
erators these expansions can be summed. Consider a 
transformation generated by 



F = A(t)I x + B(r)I y + C(t)I z 



(77) 



where r = tot and A, B and G are real functions of time. 
Let Q be an arbitrary operator that is linear in the U . A 
straight forward calculation shows that 



F, 



F,Q = a 2 F,Q 



(78) 



with 



a = Tiy? A 2 + B 2 + C 2 . 



(79) 



Therefore for any Hamiltonian linear in the spin oper- 
ators any commutation relation in (|47|l can be reduced 

to F, H or to F, F, H and the series is given by 



e ip He' lp = H - I 



F. 



+ i 
= H 



. , F,H 

2! 

- - 1 in 

F,H]--a 2 
cos a — 1 



4! a 



F,H 



F,H 



-a 

6! 



F. 



F, H 



i 

5! a 



F, 



F, H 



F,H 
sino 



o 



F, H 



(80) 



The operator F of (|77|) is linear in the spin operators and, therefore, such is also In a similar manner (|48|l can 
be summed to 



OF 
57 



dF sin a — a 



Ot 



F. 



F, 



dF 



cos a — 1 



dF 



The Hamiltonian in the new gauge is thus given by 
cos a — 1 r- r- ~ 1 1 sin a . 



G = H- 



F. 



F,H 



F,H 



. dF cos a — 1 

nuj — — h aw 5 z 

or a z 



— hio 



sin a — a 



i 7 - 



dF 



(81) 



(82) 



The problem of finding the effective Hamiltonian is thus 
reduced to finding three functions of time A(r), B(r) 
and C(t) so that l|%2l is time independent. Eq. iJS^ 
is valid for any Hamiltonian which is linear in the spin 
operators. Therefore, the problem is reduced to the so- 
lution of 3 coupled nonlinear differential equations that 
is a well defined mathematical problem. Generally this 
may be hard to do since G is not linear in terms of these 
functions. We turn now to examine the simplest case, 
that of a circularly polarized field (j75(l . 

For the spin in a circularly polarized field a perturba- 
tive solution in powers of 1/uj for F and G can be found. 
The computation is done exactly as the one of Sec. IIVI 
Thus only a brief outline of the calculation is presented. 
At the order w° 



G = H-h 



8F 1 



-luqI z + lui ( cos tI x + sinr/j, j — h— — , (83) 



dF! 



and therefore 



A = ^ (smTl x 



cos tL 



(84) 



and 



Note that here 



Go — —uqI z 



F\ , 7^ which changes some of the 

expressions obtained in Sec. IIVI 

At order w _1 a straight forward calculation leads to 



(85) 



A h dF! dF 2 lo 2 

Gi = -ncj Q — h— — I z , 

dr dr 2 



which results in 

F 2 = — ojqFi 
while 



LOqUJI 



(sin ri x 



cos tL 



cut - 
Gr = -fl, 

A similar calculation at the next order leads to 



(86) 



(87) 



(88) 



(89) 
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and to 



The Hamiltonian is then given by 



G 2 = -wqw x I x 



(90) 



The expansions for F and G are obtained by collecting 
all the terms from Eqs. H83fl - H90|) . These expansions are 
given by 



G c 
and by 



w id 2 
2lu 2 



+ OiLJ-*) I. 



(91) 



ld Q LUq 



Sid 2 



■O(oj- 



(- 



x ( sin tIt. — cos tL 



(92) 



where the subscript c denotes that this result is obtained 
for the circularly polarized field. 

An examination of Eq. (|92[1 suggests that F c may have 
the exact form 



a(u>) 
fr 



sin tI x — cos tI z 



(93) 



It turns out that F c of this form leads to a time inde- 
pendent Hamiltonian if a is chosen appropriately. Sub- 
stituting F c and H = H c of J7SJ) in JS2J| leads to 

G c = (u>i cos a — wo sin a — wsina) ^cos rl x + sinr/ y ^ 

+ (— u>o cosa — u>i sin a: + cj cosa + id) I z . (94) 

The Hamiltonian G c is time independent if 

ui i cos a — idQ sin a — w sin a = 0. (95) 

Solving for a the Hamiltonian in the new gauge, 194(1 
reduces to 



G c = [id - Jid 2 + (id +id) 2 )I z , 



(96) 



and is time independent. The quasienergies of the spin 
in a circularly polarized field are the eigenvalues of 19611 1. 
They are given by 



E x = 



(id - \Jid\ + (ido+uj) 2 ^ frs (97) 



where s = —5, —5 + 1, • • • , +5 (5 is the magnitude of the 
spin) . 

For spin S — 1/2 not only the quasienergies but also 
the quasienergy states can be computed rather easily. 
The spin operator can be represented by the Pauli ma- 
trices 



; s 1 f _ fi / -i 

^ ~ 2 I 1 ) ' l H ~ 2 I i 



t = h 

j J-z — 2 



1 

-1 



(98) 



fr _h ( -wo wie 11 
2 \ id\e ido 



(99) 



The unitary transformation U c = e~ lF " which transforms 
the eigenstates of the effective Hamiltonian (|96|l to the 
quasienergy states of (|99|l can be obtained by calculating 
the various powers of F c . For S = 1/2, 



F c =* 

2i \ e 



-e- iuJt 

iu)t Q 



(100) 



Since G c is proportional to I z its eigenstates are the 
eigenstates of I z . Thus, the quasienergy states of the 
Hamiltonian l|99|l , corresponding to quasienergies 



E± 



4<- 



arc: 



and 



u + = U c 



= U r . 



e~ luJt sin f 



(101) 



(102) 



(103) 



This is exactly the problem that was solved by Rabi [33 
and is discussed in textbooks • The physical quantity 
of interest is typically the amount of spins that flip (if 
all spins are polarized initially), rather than the Floquct 
states. We note that the term \/ui 2 + {ujq + id) 2 in the 
expression for the quasienergies is the Rabi frequency. It 
is the frequency of oscillations of these "spin flips" . 

For the spin in a linearly polarized field a perturbative 
computation of G leads to 



G, 



-id 



Aid 2 



,3, ,2 



Aid 4 



idQld\ 

6Aid 4 



0{id~ 



. M,,1} 

To the best of our knowledge an exact expression for 
the quasienergies of this system is not known. If one 
substitutes F of the form Q77JI in l)82[l. the problem of 
finding the effective time independent Hamiltonian re- 
duces to the problem of finding the three functions of 
time, A(t), B(t) and C(r), so that the new Hamiltonian 
is time independent. These satisfy first order nonlinear 
differential equations. Typically solutions to such equa- 
tions exist but it is not easy to find them explicitly. It is 
possible to choose parameters so that also the exact G c 
is proportional to I z . The approximate effective Hamil- 
tonian l|104|) can be compared to previously published 
results. While we are not aware of any 1/uj expansion 
for the quasienergies, some expansions in the strength of 
the time dependent field have been published. If one ex- 
amines, for instance, the expansion given by Eqs. (2.10) 
and App. A of 01 (which is valid for spin 5 = 1/2) and 
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expands it in powers of 1 /u> one obtains the quasienergies 

of cni. 

In this section we have studied some problems involv- 
ing spins in crossed constant and time dependent mag- 
netic fields. We have shown that the perturbation theory 
presented in Sec. IIVI can be used for such systems. For 
a circularly polarized field we were able to compute the 
quasienergies exactly, in agreement with previously pub- 
lished results. 



VII. SCATTERING FROM AN OSCILLATING 
GAUSSIAN POTENTIAL 

The systems considered in Sees. Ivl and I VII are simple, 
in the sense that the spectrum of the effective Hamilto- 
nian G is discrete and simply related to the one of the 
time independent part of the original Hamiltonian Hq. 
Moreover, for these examples also the eigenf unctions of 
these Hamiltonians are simply related. It is of interest to 
examine examples that are more complicated and where 
such simple relations cannot be found. In this section we 
examine such a system, the oscillating Gaussian, where 



an additional difference is that the spectrum is continu- 
ous and one is interested in scattering states. 

Consider a system which consists of a particle that 
interacts with an oscillatory Gaussian potential. The 
Hamiltonian is given by 

H = ^-+je-> 3x2 cos{u J t). (105) 

The system is of interest for two reasons. First, when 
x — > oo the potential vanishes and therefore one expects 
to find scattering quasienergy states. Second, the aver- 
age of the potential vanishes, namely Vq(x) = 0, conse- 
quently any interesting effect is due to the rapidly os- 
cillating potential. This system describes trapping by 
an oscillating field, a phenomenon that is of physical in- 
terest. The physical properties of this system and the 
numerical methods used to analyze it are discussed else- 
where j^jj . Here we only state briefly the results that are 
related to the properties of the effective Hamiltonian. 

The effective Hamiltonian (fH2"f . that corresponds to 
(ITOSJ, is 



G = ^ — h Prfe-*" + (1 - 2^f e^f - ±^X^ - l) (3 - 2^) e^p 

- H (-9 + 99/3x 2 - 114/3 V + U(3 3 x 6 ) e^ &3? + 0(i^ 5 ). (106) 



We examine separately the leading correction due to the 
oscillating field, that is given by the Hamiltonian 

G< 2 > = G + G x + G 2 = f- + ^J^ e -^ 2 , (i 07 ) 
2m mw z 

where the error is of order w -4 . It has a simple physical 
meaning as the Hamiltonian of a double barrier potential 
and its spectrum is continuous. 

Since the effective potential of Eq. (|107|) is a double 
barrier one expects to find that this system exhibits res- 
onances. These resonances describe long lived unstable 
states. Each resonance is characterized by a complex en- 
ergy E — iT/2. The real part E is the location of the 
resonance while Y is the width which is inversely propor- 
tional to the lifetime. For a review on relevant proper- 
ties of resonances and useful methods to compute them 
see H2. 

For any resonance of (|106f> and l|107|) it is natural to 
look for the corresponding resonance of the time depen- 
dent original Hamiltonian 1|105|) . More precisely, one 
looks for the resonances of the Floquet Hamiltonian Ij33(l 
with H of l|105|) . This is done numerically using a com- 
bination of the (t, t') method and complex scaling |4^ . 

The energy Eq and the width Tq of the lowest (cor- 



responding to the smallest real part Eq) quasienergy- 
resonance of (I105|l are compared with the lowest reso- 
nance of the effective Hamiltonians (I106f) and 11071) in 
Figs. Hand [21 It is clear that for large frequencies there 
is excellent agreement between the resonance of the time 
dependent Hamiltonian i|105|) and the ones of the effec- 
tive Hamiltonians (|106|) and (|107|) . At low frequencies 
the location and width of the exact resonance differ from 
those of the effective Hamiltonian. The deviation for the 
order lu~ a Hamiltonian (|106fl is large indicating that the 
expansion is asymptotic. This is expected since the per- 
turbation theory developed in Sec. IIVI assumes high fre- 
quencies. A more complete study of this specific system 
and a discussion regarding the physical implications of 
this resonance are given in |4l| . 

In this section we demonstrated that the effective 
Hamiltonian G can be used to obtain some physical prop- 
erties of systems that are more complicated than those 
presented in Sees. [V] and IVII In particular, the reso- 
nances of a periodic time dependent system were found 
to be given by the resonances of the corresponding time 
independent effective Hamiltonian. Resonances for os- 
cillating barriers were computed numerically in |l7l Il8| . 
The calculation of the present section demonstrates their 
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Time dependent 
Effective - to 2 
Effective - CD 




FIG. 1: The energy Eo of the lowest quasienergy resonance of 
the Hamiltonian 110511 as a function of the driving frequency 
(solid line) , compared to the lowest resonance of the effective 
Hamiltonians 110711 -dashed line, and 110611 -dotted line, for 
7 = 9 and j3 — 0.02. The "atomic units" m = Ti = e = 1 are 
used here. 
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FIG. 2: Same as Fig. for To, the width of the lowest reso- 
nance. 



physical origin. 



VIII. SUMMARY AND DISCUSSION 

In this paper we investigated classical and quantum 
motion in high frequency fields. The classical motion 
can be treated by separation of time scales. In Sec. |H] 
this motion is separated into a slow part and a fast part, 
which consists of rapid oscillations around the slow part. 
The fast part and the resulting equation for the slow 
motion are solved perturbatively to order oj~ . This per- 
turbation series is a generalization of the calculation pre- 
sented in Mechanics by Landau and Lifshitz ■ We also 
demonstrated that this perturbation theory is equivalent 
(to order w~ 4 ) to the standard mathematical method of 
multiple time scales analysis |43j . which is more com- 
plicated. The resulting equation for the slow motion is 



found to result from a time independent Hamiltonian. 

Following a review of Floquet theory in Sec. lIIII the cor- 
responding effective quantum Hamiltonian is computed 
explicitly, using a high frequency perturbation theory up 
to order u> , in Sec. II VI The resulting Hamiltonian l|62() 
is rather simple. Its classical limit is the classical effective 
Hamiltonian (|3U|) . This effective Hamiltonian is there- 
fore a generalization of the classical results of Kapitza 
and Landau and Lifshitz 7] to quantum mechanics. In 

the classical limit the unitary gauge transformation e~ lF 
reduces to a canonical transformation generated by the 
classical limit of — h~F (see (|38H and 13]). The Hamilto- 
nian (|3U|I for the slow variables was obtained in a natural 
way in Sec. [H] The freedom in the choice of gauge was 
used to choose F so that in the classical limit the effec- 
tive quantum Hamiltonian (|62[) reduces to H3U|) . Con- 
sequently the classical limit of — fiF is the generating 
function of the canonical transformation from the orig- 
inal time dependent Hamiltonian to the time indepen- 
dent Hamiltonian l|3U|) . This limit explains the fact that 
the classical dynamics of the slow coordinate X is gener- 
ated by a Hamiltonian. Using the freedom in the choice 
of gauge one can generate Hamiltonians that differ from 
(|3U[) and (|62[l but are related to them by canonical and 
gauge transformations. The present work extends [L3I ] 
to general driving potentials and is not restricted to the 
driving (6a) of [l3| . The perturbation theory which was 
developed can, in principle, be used to compute it to any 
given order in 1/lo. This is a significant extension be- 
yond |li| in the spirit of separation of time scales 
that enables a systematic expansion in powers of u) . 
For this, the requirement (23) of 01 is avoided and the 
expansion can be performed for any driving potential V± 
and is not restricted to driving forces that are uniform 
in space. It should be emphasized that this perturba- 
tion theory is an expansion in 1 /u> and not in powers of 
the time dependent potential. The potential V\ does not 
have to be small in order to obtain a good approximation 
of the original system. 

Several examples were discussed. The driven harmonic 
oscillator and the spin in a rotating magnetic field are 
simple, exactly solvable examples, that were used as 
demonstrations for the method. For another system, the 
oscillating Gaussian, we showed numerically that its low- 
est resonance is given by the resonance of the correspond- 
ing effective Hamiltonian I|ri2l) of Scc. lVIII Thus, for time 
dependent traps, such as the atomic billiards discussed 
earlier, the time independent effective Hamiltonian can 
be used to compute the resonances and the lifetimes of 
particles in these traps. 

While the examples presented in this work indicate 
that this effective potential is a meaningful concept and 
is also useful for calculations, there are points that re- 
quire further research. The convergence properties of 
the expansions for F and G are not clear. There 
may be situations in which the perturbation theory fails 
to converge or where smaller than any power corrections 
are of physical importance. Consider for example a sys- 
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tern similar to the oscillating Gaussian where V\ vanishes 
outside some finite domain. In this case one may expect 
to find scattering states as the eigenstates of G. In the 
limit of high energies such a state may be roughly ap- 
proximated by the plane wave e lpx / h (with large p) . The 
leading momentum contributions to F are of the form 



1 



[Vi. 



1 



hmuj 2 



(2)t 



[v(]p- 



1 



(3)t 



[V{']f 



hm 2 uj 3 



Twm?io A 



(3) 



P 



as one can see from (|66|) with (|59|) and (|C7|) . While the 
general structure of the series for F is unknown, it seems 
that for states with momentum which scales as to a with 
a > 1 the series do not converge. 

Even taking into account all the questions concerning 
the validity and convergence of the perturbation expan- 
sions presented in this work we find the effective potential 
to be useful for the following reasons. The perturbation 
theory leads to a time independent effective Hamiltonian. 
Physicists, who are used to work with time independent 
systems, have developed an intuition for such systems, 
and thus the effective Hamiltonian may give physical in- 
sight that is absent when examining the corresponding 
time dependent problem. In addition, all the calculations 
used to obtain the effective Hamiltonian are straight for- 
ward. There are no differential equations or any compli- 
cated iterative schemes that may appear in more sophis- 
ticated perturbation theories. For comparison see |44| . 
Finally, one may use all the well developed techniques 
for time independent quantum systems to compute the 
eigenvalues of G. In particular, one can use time inde- 
pendent perturbation theory in the case where the eigen- 
values and eigenstates of Go are known. 

The effective Hamiltonian l|62|) can be useful to predict 
the qualitative behavior without complicated numerical 
calculations. Assume first that the frequency is suffi- 
ciently high so that only terms of order up to 1/ui 2 should 

be included and denote V2(x) = ^tu? (J T ^l') • Let ^° 
take the form depicted in Fig.^Ji. It exhibits resonances. 
A natural question is whether the line width increases or 
decreases as a result of the time dependent potential. It 
is clear that for a situation of Fig. |3t> the line width de- 
creases, since the particle has to tunnel through effective 
barriers which are higher than those of Vb because the 
effective potential is V e f / = Vb + V2 , while for Fig. ^ 
the line width (typically) increases since the energy of 
the resonances is shifted upwards by the time dependent 
perturbation. Numerical calculations of the type pre- 
sented in Sec. IVHI and in |l7j and [l8j should confirm 
these results. If the terms of order u> are important 
more subtle considerations concerning the kinetic energy 
(p 2 ) /2m and the semiclassical limit are required. 

Since there are some physical situations where one ex- 
pects that the perturbation theory fails it is of impor- 
tance to find ways to improve it. The goal is to be able 
to describe terms that are smaller than any power in 




FIG. 3: Qualitative effect of driving on trapping times: a) 
the average potential Vb, b) the position dependence of the 
driving potential Vi (dashed line) results in the effective po- 
tential V e ff (solid line), c) as in b) with driving mainly inside 
the trap. The units are arbitrary. 



the perturbations. This might be possible using super- 
convergent perturbation theory that was recently applied 
to time dependent quantum systems |45| . In this pertur- 
bation theory the small parameter is the size of the time 
dependent potential. It is of interest to modify it to a 
perturbation theory in 1/lo. Note that while this super- 
convergent perturbation schemes have superior conver- 
gence properties, they may turn out to be complicated 
for explicit calculations. Thus, one may lose the main 
advantage of the perturbation theory presented in this 
work, its simplicity. In addition, it is of interest to gener- 
alize the results of this work to systems of higher dimen- 
sion. Such a generalization should be straight forward. 

In conclusion, we have investigated the dynamics of 
high frequency driven classical and quantum systems. 
High frequency perturbation theory was used to obtain 
an effective time independent Hamiltonian for the slow 
part of the classical and quantum motion. For quan- 
tum systems, the spectrum of this Hamiltonian is the 
quasienergy spectrum of the time dependent system. 
This effective Hamiltonian is computed in a high fre- 
quency systematic perturbation theory. It is demon- 
strated that the effective Hamiltonian gives the exact 
quasienergies and quasienergy states of some simple ex- 
amples as well as the lowest resonance (including the life- 
time) for a time dependent atom trap. 
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APPENDIX A: 



MULTIPLE TIME SCALES 
ANALYSIS 



The derivation given in Sec. HP was, in some sense, not 
explicitly consistent. For example the slow motion was 
not solved by expanding its coordinates in orders of 1/lo. 
Terms of order l/ui n included all contributions up to that 
order, and also some contributions of higher order. It is 
of interest to show that the same result can be obtained 
using a standard method, namely the method of multiple 
time scales analysis 0. In this appendix we show how 
to derive the equations of motion of the slow dynamics 
using multiple time scales analysis. We will present only 
the first few orders in l/u> since this method turns out 
to be more complicated than the method used in Sec.lTTl 
While the two methods differ in details they are equiva- 
lent and lead to the same results, when consistently ex- 
panded order by order in l/u>. 

In order to use multiple time scales analysis it is conve- 
nient to transform equation to a standard form [43| . 
This can be done by defining r = cut, e = — and 
y = jjj — uj^. The first order equations of motion are 
then given by 



-j^Z= ef(r,z) 



where 



and 



/ = 



-±W(x)+V{(x,t)) 



(Al) 



(A2) 



(A3) 



Then one introduces the following expansion of the solu- 
tion 



(A4) 



n=0 



where t = er and r are treated as independent variables. 
Using this expansion together with 



d _ d 
dr dr 



d_ 



results in 



dz dz 7( _> _ 

+ e— + e— H = ef (r, z + e%\ 

or at ot 



(A5) 



(A6) 



The solution is obtained by expanding /, matching pow- 
ers of e, and solving for z order by order. This is the 



standard multiple time scales analysis, see |43| for a de- 
tailed description of the method. (Note that in the 
role of t and r is opposite to the one in the present work) . 
We proceed to solve the first few orders in e. 



At order e° the leading order of l|A3|l results in 



dx 
dr 

dyo 
dr 



(A7) 



Therefore xo and yo can be any functions of the slow time 
t 



x (t) 

yo(t). 



(A8) 



Note that at the leading order the solution is found to 
depend only on the slow time scale, as expected. We will 
denote the r independent slow part of the solution by Xi 
and §i at any order. Additional conditions ijAlOfl on xq 
and j/o w iH be obtained from the requirement that the 
solution is not secular at the next order. 
At order e (| A3|> results in 



dx _ , , 

-^-L V {(x )-±-V{(x ,t), (A9) 
at m m 



dxi 

dm 

dr 



where we have used (|A8ll . We are interested in solutions 
of l|A9|l that are not secular, namely, that do not grow 
with the fast time scale r. The RHS of equation (|A9(I can 
be decomposed into periodic functions of r (with vanish- 
ing average) and r independent terms. Any r indepen- 
dent term will result in a secular contribution. Therefore, 
to avoid such terms we demand 



dx 
dt 

dy 
dt 



--Vi{xo{t)). 



(A10) 



This is the leading order of the equation that governs the 
slow time scale. The non secular equation of order e is 
now 

IT =° 

OT 

dyi 
dr 

and its solution is 



-Vl(x ,r) 



(AH) 



xi(r,t) = xi(t) 

Vi(r,t) = -- [ [V;(xo(t),T)}+ m (t). (A12) 

TO J 

This process can be repeated order after order. At 
each order all the terms with the same power of e are 
gathered. The resulting equation may be secular and 
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therefore an additional condition on the slow part of the 
solution is enforced. Then, one can solve for x and y at 
that order. The calculation is rather tedious and will not 
be presented here. We only give the secularity conditions 
that result when the next two orders are computed: 



and 



dx i 

~dt 

dyi 

at 



yi(t) 
1 



siWvtfW*)), 



dx 2 

dm 
dt 



(A13) 



(A14) 



-Lv['(x (t),T) J T [V{(xo(t),r)} 

1/(3), 



-x 2 (t)Vj'(x (t)) 
m 



2m 



Equations (|AT0)l . (|A13)) and (|Al4|l are the first orders 
of an expansion of the following system of equations 



dX 
~dt 
dY 
~dt 



Y(t) 
1 



m 



—V{'(X(t),r) 



(2)7 



[V{(X(t),r)] (A15) 



in powers of e where 



X(t) 
Y(t) 



x (t) 



m{t)- 



%(t) - 



(A16) 



Equations (|A15I) are accurate to order e 2 and are identi- 
cal to the slow equation (|22|l . We have used the multiple 
time scales analysis up to order e 4 and found that the re- 
sulting slow equations analogous to i|A15(l are equivalent 
to equation (|27|l in Sec.[n] with the identification X = X 
and Y = X. 

The method presented in Sec. [H] and multiple time 
scales analysis used in this appendix have a similar struc- 
ture. Both separate the motion into r dependent and in- 
dependent parts and both lead to equations for the slow 
motion that have to be satisfied to avoid solutions that 
grow with t. The main difference between the methods is 
that when using multiple time scales the slow coordinate 
is expanded in powers of e = 1/ui. This leads to a large 
number of terms that result of the fact that the functions 
in the slow equation are evaluated at Sio(t) rather than 
at X(t). Consequently the derivation in Sec. |n]is much 
simpler than the one presented in this appendix. This 
is also the reason that here we did not present explicitly 
the calculation using multiple time scales analysis up to 
order e 4 . 



APPENDIX B: SOME USEFUL OPERATOR 
RELATIONS 

In this appendix some known relations (see for exam- 
ple ^(|) involving operators and exponentials of other 
operators will be presented for completeness. Let us de- 
fine 

C{x) = e xA Be~ xA . (Bl) 
It is clear that C(0) = B. Differentiation shows that 



dC 
dx 



AC 



(B2) 



Therefore C(x) has the following expansion in powers of 
x 



C 



B + x 
1 



A 



A, B 
A, 



1 

+ 2 : 
A,B 



A. 



A,B 



(B3) 



Substitution of x — 1, A = iF and B = H results in l|T7Jl . 
One can use the expansion Ij47(l to find an expansion 

~ lF . First, since e lF e~ lF = 1, 




-iF 



de 



-iF 



= -e iF —e- iP + —, 
dr dr ' 



(B4) 

where in the last expression the derivative ^ in front of 

e~ zF operates not only on the exponential but also on 
any operator or function that appears on its right. Using 

(|B3)l with B = £ and A = iF in (JB4J) yields 



de 1 



-,-iF 



f A" 


1 




' f A" 




'dr 




F, 


'dr 




+ 2 





I 

3! 



F. 



F- 



Ot 



(B5) 



Substitution of 



0_ 



dF 



(B6) 



in HB5[1 results in |@SJ|. 

If the operators A and B commute with their commu- 



tator 



A,B 



then m 



e A+B = e-il A ' B le A e*. 



(B7) 



This is a special case of the Campbell-Baker-Hausdorff 
formula [47| . 
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APPENDIX C: SOME HIGH ORDERS OF THE of order w 3 found from gUl gHJ is 
EFFECTIVE HAMILTONIAN 



In this appendix we will outline the steps that should 
be taken in order to compute G to order uj~ a . The term 



G, 



F 3 ,H 



F 2 



F U H 



ih 

h 

6 



F 2 , 



d_F2 
dr 



~2 



Fa, 



8Fi 



F u 
h r 

6 



F 2 ,H 



F\ 



Fu 



F\,H 



- h 



8F 4 ih 

"d7 ~ IT 



Fi, 



0F3 
dr 



Fi 





fi 9Fx 




ih 












Fi, 




Fl '^7 




+ 24 



Fi, 



Fi 



F 



dFi 



dF 2 
dr 



h 
6 



Fx, 



F 2 



dFi 



(CI) 



It can be further simplified by taking into considera- 
tion that some of the commutators vanish, for instance, 



Fi, 



Fo 



dFi 
dr 



= since there is only one derivative 
(with respect to x) and two commutations. Equations 
and (|58Jl for F% and F3 can be used to eliminate 
the commutation relations that involve H. After some 
manipulations (|C1|) reduces to 



G, 



F3, Gq 



dFi 

OT 



ih 

~2 



Fo 



8F\ 
dr 



(C2) 



F4 will be chosen so that the periodic terms (with van- 
ishing average) in (|C2(I are eliminated. To examine the 
contributions to G3 which remain, one can just average 
J2I) over a period. This leads to 




(2)r . 



(3) 



(C3) 



The freedom to choose / 3 can be used to cause G3 to 
vanish. We choose /a to be a function of x so that 



1 f T r (2)r 
Using ljC4() in (|C3|) indeed results in 



(C4) 



G 3 = 0. (C5) 
Thus we will use this choice of gauge. With this choice of 
f 3 (and therefore of F3), we require that F4 will satisfy 



dF A _ i 
~dr ~ h 



F3, Gc 



dFi 



(C6) 



Equation (|C6|I can be integrated in order to compute F4. 
Using JSTJ, iJSSJl and QCty leads to 



fi = -t± [ i4)T \ V m «L _ 3^1 [ i4)T \v^] *- - ™L1 [ (A)T M 5 >1 ± - f (4)T M 6 >1 (C7) 
4 m 3 J I 1 J fe 3 2m 3 y L 1 J dx 2 4m 3 J L 1 J dx 8m 3 7 L 1 J K ' 

- 5? (/' w + /' + 2 f w s + 2 /' w /** 1*1 + 2 f K1 5) + ><<*•»• 

where fu is a hermitian time independent operator (that depends only on x and p) that will be determined at the 
next stage while 

H3)t 7(3)7 

^(x,r) = VI J [V{']-V(J [V{>] 

A3)r ,(3)t 

V 3 (x,r) = V "j [V(]+3V^J [V"} . (C8) 
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The operator fa plays here a role similar to the one of fa. Note that V\ = V 2 = V3 = 0. 

We turn to compute the next order of G, order w~ 4 . This is the last order that will be considered explicitly in this 
work. The terms of order w~ 4 in (|44Jl - 1)48(1 are given by 



G 4 = 



where 



h 

120 



4 



1- i fi 1 [ 
~2 £l "6 £2+ 24 



Fx 



Fi, 



* dF 5 i?l KA 

~dV ~ ~2 1 



-M 2 



24 



M 3 







Fx, 


Fx, 



Mi = 



Mi 



Fi, 
Fi, 

A, 

A, 
A, 



Fo,H 



Fn 



Fx, 

dFx 
dr 

Fx 



F 2 ,H 



F 2 ,H 
A 



F 



3, 



dF 3 
' dr 



- dF 2 



8F2 
dr 

Fx 



F 2 , 



F 2 , 

A, 



F 3 
Fx,H 

dr 



Fi 



F 2 



dFx 
~d7 

dFx 
~d7 



Fx,H 



Fx, 



F, 



F 2 , 
dP 4 



Fx, 



Fx,H 



(C9) 

(CIO) 
(Cll) 

(C12) 



Fx 



dFx 



Fx, 



dF 2 
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and 
M 3 = 







p dF 2 ' 












a dFx' 












A dFx' 












a dFx' 






Fx, 


Fx, 


Fu ^7 






+ 


Fx, 


Fx, 








+ 


Fx, 


A, 








+ 


A, 


Fx, 









(C13) 



(C14) 



Many commutation relations vanish since F 2 contains one derivative (see (|54|l 'l and Fx contains no derivatives. Re- 
peated application of the commutation relations results in the vanishing of such relations. The commutation relations 



of H (except 



F 4 ,H 



can be simplified to 



G& 



Fi, Go 



can be eliminated using (|50|l . (|53[l . (|58() and (|C6|I . After a tedious calculation, equation (|C9|) 
t dP 5 ih 



F 



Fo 



dFx 



dr 



Fa 



dFx 



h 
12 



Fx, 



Fo 



dFx 
~d7 



n 
12 



F ? 



F ? 



dFx 



(C15) 



Our goal is to obtain an explicit expression for G4 as was done for lower orders. Note that we do not have to compute 
P5 explicitly for this purpose since P5 is chosen in such a way that it will cancel the time dependent part of IJC15I) and 
it is required only for the calculation of higher order terms of G. The terms that are not canceled are just the time 
independent terms on the RHS of (|C15|> . They are obtained simply by averaging (|C15J| over a period resulting in 



Ga 



fa, Gq 



ifr 
~2 



F 4 , 



dFx 
dr 



12 



Fx, 



F 3 , 



dFx 
dr 



12 



F 2 , 



F 2 , 



dFx 
dr 



(C16) 



Substitution of the expressions for Fi and Go leads after some straight forward but tedious calculations to 

-2 

G 4 




1 f / /■(2)r \ 2 r (2)r r (2)r ,r ,(3)r ) 

3^ V " [ J [V(] ) + V{ J lV(] J [V{ ' ] ~ 2V( J [V(] J [V " ] 



(C17) 
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Using integration by parts it is possible to show that 



Vi 



,(2)r ,(2)r r r ,(3)t i / r (2)r \ z 

i y y - 2v? y [Vj] y [vc] = -v" f y ^ j . 

This will simplify the last term of (|C17|I . The operator fi is not determined yet. We choose it of the form: 

h = 9(x)p + pg(x) = 2-g{x)^- + -g'{x) 
i ox i 

where g(x) is a function that should be determined. In terms of the function g(x) we find 

-2 



G 4 = 2hgV( 




1 / /-^ T \ /- (2)r /- (2)T 

a^My raj +3^y vi] j [ V {>] 



(C18) 



(C19) 



(C20) 



We choose g(x) so that in the classical limit the effective Hamiltonian G reduces to its classical counterpart <|30ll . It 
takes the form 



Q r(2)T A2)r 



(C21) 



The resulting operator G4 is obtained by substituting IJC21I) in (|C20(I . It results in the term of order u> 4 in the 
effective Hamiltonian (|fi2l) : 



1 / r Wr 

G 4 = ^U " / [V{] 



2m 



1 ( /-(2)r V t 2 



2m 3 



(2)t 



^1 

9 a; 2 



(2)r 

6 / [V/' 



(2)7 



(3) 



d_ 3 /- (2) 

+ 2 



(2)r 
/ 



(2)r r 



(3) 



(C22) 



This is the highest order of G that is computed here. 
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